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We present an efficient graphene-based photodetector with two Fabri-Perot cavities. It is shown that the 
absorption can reach almost 100% around a given frequency, which is determined by the two-cavity lengths. 
It is also shown that hysteresis in the absorbance is possible, with the transmittance amplitude of the mirrors 
working as an external driving field. The role of non-linear contributions to the optical susceptibility of graphene 
is discussed. 
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I. INTRODUCTION 

Exploring the optical properties of graphene^] for photode- 
tection is one of the most promising applications of graphene. 4 
Graphene has no gap and its conductivity is essentially inde- 
pendent of frequency 5 9 for photon energies up to 2 eV. These 
properties, combined with the intrinsic chemical and mechan- 
ical stability of graphene, pave the way for broad band opto- 
electronics. 

Depending on the authors, graphene is characterized as pre- 
senting remarkable high absorption^ or weak absorptionpl 
Indisputable however is the fact that pristine graphene ab- 
sorbs about 2.3% of the light impinging on it. A consider- 
able frequency dependence of the absorption appears as the 
photon energy approaches ~ 4.7 eV, due to the combined 
effects of the Van Hove singularity in graphene 's electronic 
7T— spectrum and excitonic many-body effects P^l The value of 
2.3% can also be interpreted as the probability for photon ab- 
sorption in a single passing through the material. 

To enhance the absorption of graphene several mech- 
anisms have been proposed, ranging from hybrid mate- 
rials, containing carbon, boron, and nitrogenp-^ n ano- 
patterning of a graphene sheeP^, to strain engineering! \UiU\ 
and plasmonics! 15 * 16 ! In the latter case, micro-sized ribbons 
patterned in a single graphene sheetP^ and metallic arrays on 
top of graphene, 18 lead to an enhancement of the near field, 
thus increasing light absorption and producing larger photo- 
currents compared to the case of pristine graphene. Current 
plasmonics-based approaches are limited to specific spectral 
bands. In what concerns hybrid materials for photonic appli- 
cations, there is still a long way to go before these become 
possibleP^^ 

Photodetection, depending on the type of application, may 
require efficient absorption of light in a narrow spectral band. 
Then, it is conceivable to explore Fabry-Perot interference for 
producing an efficient photodetector tailored for a specific ap- 
plication. The concept is simple: in a Fabry-Perot interfer- 
ometer a photon may be trapped inside the cavity for a long 
time, undergoing many round trips before leaving it; indeed, 
if a material able to absorb photons is introduced inside the 
optical cavity, without significantly changing the cavity's fi- 



nesse, as is the case of graphene, most of the photons of the 
right frequency entering the optical cavity will be absorbed. 

For a combined cavity-graphene system, it is important to 
quantity the magnitude of non-linear optical effects, ensur- 
ing whether linear response theory can be used for describ- 
ing the absorption process inside the cavity. Thus, in the 
present work, we discuss the non-linear optical susceptibility 
of graphene first, presenting later the single and double cavity 
photodetector s. 

This article is organized as follows: in Sec. [II] we introduce 
B loch's equations for graphene and compute both the linear 
and non-linear optical susceptibility. As aforementioned, the 
calculation of the non-linear part is essential for a critical anal- 
ysis of the optical response of graphene. In Sec. [Ill] we give 



the power series solution of B loch's equation and discuss the 



validity of perturbation theory. In Sec. IV we introduce the 
mathematical description of a graphene-based photodetector. 
Having shown that non-linear contributions to the optical sus- 
ceptibility are relevant only at very high field intensities, we 
employ linear response theory to describe the graphene-based 
photodetector with two coupled optical cavities. 



II. DERIVATION OF BLOCH EQUATIONS OF MOTION 

The calculation of the optical properties of a given material 
can be obtained from the solution of B loch's equations E3E3 
Below, we derive B loch's differential equations for graphene 
and give its solution for an incoming electromagnetic plane 
wave impinging on the material. 

The Hamiltonian of the electrons in graphene, in the pres- 
ence of an electromagnetic field, reads (spin index implicit): 

H = ^2 E c(k)a\ M a Cik + E v {k)a\ M a Vik 

k 

+ v F eA(t) ^2(d x cv ^ k al k a Vjk + d x vc k a\ k a cM ) , (1) 

k 

where E c / V (k) = ±vphk, vf = 3tao/his the Fermi velocity 
(t ~ 2.7 eV and ao = 1.4 A are the hoping integral and the 



2 



carbon-carbon distance, respectively), a£ , fe is the creation 
operator of an electron with wave number k in the conduc- 
tion/valence band, 



d x 



(2) 



is the matrix element of the dipole operator {d x cv k = ism 6), 
A(t) = A sin(u;t) is the vector potential for a linearly po- 
larized electromagnetic plane wave, e > is the elementary 
charge, u is the frequency of light, and 



Since we have E(t) 
write A(t) as 



— arctan — ^ . 

dA(t) 
~ dt 



v ' 2uj 2iuj 



(3) 

wAq cos(ut) , we can 
(4) 



where Eq is the intensity of the electric field. 

At the heart of the present approach is Heisenberg's equa- 
tion of motion for the polarization operator P vc ,k = <4 fe a c,fe' 
namely, 



(5) 



The explicit form of the equation of motion is 

- ^l2,kPvc,k = -[E c (k) - E v (k)]P VC) k 

+ v F eA(t)(a ] c k a cM - aJ jfc a Vjfc )d^ jfc . (6) 

In the above, 72, & is the phenomenological relaxation rate of 
the polarization. In addition, we need the equation of motion 
for the number operator, h\^, both in the conduction and va- 
lence bands, reading 

- ~ « ft 7i,fc^c,fc = ev F A(t)D VC)k , (7) 

.^dn^k _ ifr lx h fo v h = ev F A(t)D* vc k , (8) 



where 



dt 



7~) JCC p J# p 

J-^VC.k — U J vc ± vc k U J cv ± cv k ; 



(9) 



and 71^ is the phenomenological relaxation rate for the oc- 
cupation number of a state k. It is convenient to define the 
population operator, Nk = n c ^k — n v ^k, whose average obeys 
the following differential equation 

H Jt^ k + Hll ^ k = 2eVF 0A (t)(Pvc,k + Pcv,k) • (10) 

We denote the average of Nk, Pvc,k, and P cv ,k, by Nk, P V c,k, 
and P cv ,k, respectively. The latter averages obey the follow- 
ing set of linear, first order, differential equations (known as 
B loch's equations): 

d Pvc,k + 72,kP vc ,k = -itkPvck ~ d k A(t)N k , (11) 



dt 
d 

It 



Pcv,k + l2M P cvM = iekP C v,k ~ d k A(i)N k , (12) 



d 
dt 



N k +~fi,kNk = 2d k A(t)(P v 

c,k H~ Pcv,k , (13) 



where dk — v F esmO/h and e*, = 2v F k. We note that the 
differential equation for P cv ,k is redundant, since P cv ,k = 
[Pvc,k]^ • The differential equations are solved together with 
the initial conditions (t = —00.): P vc ,k = and Nk = 
no = f[E c (k)] — f[E v (k)], where f(x) is the equilibrium 
Fermi distribution. The condition no = — 1 applies to neutral 
graphene at zero temperature. 

We note in passing that in the absence of relaxation mech- 
anisms, i.e., 7i, fc = 72, k = 0, the quantity N% + 4:P VCjk P C v,k 
is a constant of motion. We also note that there is no funda- 
mental reason why 71^ should be equal to 72, albeit they 
generally are of the same order of magnitude; in order to 
simplify the mathematical expressions, we assume below that 
7i,fe = 72, k = 7- This procedure is justified since making 
7i, fe 7^ 72, fe does not alter the final qualitative conclusions. 



III. POWER SERIES SOLUTION TO BLOCH'S 
EQUATIONS 

In the present section we obtain the solution of B loch's 
equations. It is convenient to employ the shorthand nota- 
tion, P VCyk X(t) = X, N k = N(t) = N, c k = e, and 
dkA(t) — a(t). Using this notation, B loch's equations have 
the form 



d_ 

dt 
d_ 

dt 



X + 7 X = -ieX - a(t)N , 



(14) 



N + 7X = 2a(t)(X + X*) = 4a(£)J?X . (15) 



We now assume a power series solution for X and 



a(t) —> Xa(t) , 

X(t) -> \X!(t) 

N(t) -> n 



\ 2 x 2 (t) + \ 3 x 3 (t) + ... 
Am(t) + A 2 n 2 (t) + A 3 n 3 (t)- 



(16) 
(17) 
... , (18) 



where A is a bookkeeping of the power of the amplitude of the 
electric field; it is useful to use the notation n& = no. Intro- 
ducing the series expansions ( [T6| ), ( [T7] ), and ([18]) in Eqs. ([14]) 
and ( [T5] ) it is simple to see that only odd powers, X2m+u of 
the polarization are non-zero, whereas for the population only 
even powers, n2 m , are finite, with m an integer number, in- 
cluding zero. 

Two relevant dimensionless parameters involving the inten- 
sity of the incoming field, W%, are: 



/3 7 



27 Wjalt 2 
4 H 3 u 2 j 2 



and 



27 



malt 2 



4 /i 3 w 2 (w 2 +7 2 /2) 



27 Wialt 2 

j 

4 ft 3 u; 4 



(19) 



(20) 



where Wi — EqCqc/2, a = e 2 /(Aneohc) is the fine structure 
constant, and we also have assumed u ^> 7 in When ei- 
ther /3 7 > 1 or 13^ > 1, the perturbative solution breaks down 
and the full series has to be resumed. The choice of prefactors 
in /3 7 and will be apparent later in the text. Numerically, 
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the intensities setting the limit of validity of perturbation the 
ory are 

,GW 



w il7 = (M (fry) x io d 



from /3 7 =1, and 



= (M 4 x 10 ; 



GW 



cm' 



(21) 



(22) 



from /3a;=l, with fouj and fry expressed in electron-volt; for 
graphene we have fry ~ 10 meV. Taking a representative 
value of fro> ~ 0.5 eV, we obtain 



W i7 ^2.5xl0- 2 ^ 



W i<w - 60- 



GW 
^2 • 



(23) 
(24) 



It should be noted the three orders of magnitude difference 
between the two cases. 



A. Linear optical susceptibility 

The calculation of the optical susceptibility of graphene can 
be made for any frequency value P^32U Q n fa e ther hand, 
Hamiltonian ([T]) is valid up to energies of the order of 1 eV, 
which translates into photon frequencies of the order of 2 eV. 
Hence, both for illustrating the method and describing how the 
photodetector works, the Dirac cone approximation suffices 
for our purposes. 

The solution for the linear polarization, x\(t), is obtained 
from 

±1 + (ie + 7)^1 = -a(t)n , (25) 

which is easily solved by the integrating factor e^ e+7 ^, lead- 
ing to 

Xl {t) = - e -( <c +i0* C e (* e +7K a (t f )n dt f . (26) 
j — 00 

In the particular case where A(t) is described by a sinusoidal 
function we obtain 



xi(t) = -n k d k 



En 



-iujt 



' 2lj \lj + e — 27 co> — e + 27 
In general, the total polarization is computed from 

P x = -ev F g s g v ^(P vc ,fc< c + P CVyk d* v ) 

= AP^i + A 3 P X5 3 + . . . , 



. (27) 



(28) 



where g s and g v are the spin and valley degeneracy, respec- 
tively. Recalling that P VCjk = X and P cv , k — X*, we have, 
to first order in the electric field amplitude (Eq = Aq/ilj) 



P x ,i = iev F g s g v ^ sin 0[xi (t) - x*(t)] 

k 



(29) 



with Xi( — u ) = Xii^). Considering the case of neutral 
graphene at zero temperature, we obtain for the real part of the 
optical susceptibility, 5Rxi(^) = xi» me well known value 



(30) 



dubbed the universal conductivity of graphene. 2 Given 
Eq. ( [30] ), the imaginary part of the optical susceptibility reads 
^Xi(^) = Xi — 0> as follows from the Kramers-Kronig rela- 
tion: 



1 f c 



Xi(a) 



dx . 



(31) 



Equivalent relations hold for non-linear response functions as 
well.™ 



B. Non-linear optical susceptibility 

To go beyond linear response, we have to compute how the 
population changes relatively to its initial value when the field 
is turned on. This amounts to compute ri2(t). The latter can 
be obtained from the solution of x\(t) according to 

712 (t) = 2e~^ f dt'a(t')MO + x\(t')] . (32) 



The occupancy second-order correction, 77,2 (t), is a sum of 
two contributions: 77,2 (t) = ri2 a + ^26 (t), where the first one 
is independent of time. Explicitly we have 



n 2a 



and 



, E n 



1 



1 



~ nfcd * J V(^-^) 2 +7 2 ' (a; + e fc ) 2 + 7 2 



(33) 



n 2 b(t) 



2 h u^ 



-2iout 



2uj — ( uj — e k — 



4uj 2 + 7 2 \ (lj - e k ) 2 + 7 2 



uj + e k - Z7 
(u; + e fc ) 2 +7 2 



c. c. 



(34) 



We note that ri2 a is a positive number, thus reducing the value 
of N(t) when the system is driven away from equilibrium. 
Similarly, the calculation of x$(t) follows from 

x 3 (t) = -e-^^ f e^^' a(t')[n 2a + n 2h (t')]dt' • 
j — 00 

(35) 

The quantity xs(t) is a sum two different terms x 3 (t) = 

X3a(t) +X 3b (t): 



iuit 



E 3 e 

x 3a (t) = in k d 3 k —^A(e kl uj)- — — ■ — 

4ur iuj + ie k + 7 



+ (cj — » — cj) , 
(36) 

where (lj — » — a;) is obtained from the given explicit term 
upon the replacement uj — )> — cj, and 



^3 gicjt 

^36<» ~ -in k d k — \B(e k ,uj)- — — — 

4ur zcj + zefc + 7 



+ (lj -> -uj) . 

(37) 
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Terms proportional to e ±l3u)t corres pond to three photon ab- processes, as represented in Fig. [T] Each of these processes 

sorption and were neglected in Eq. ([37]) PS We have also de- excite different electronic states. 

fined For sake of completeness, we give the formula for the op- 

2 2 ^ical susceptibility due to three photon-absorption processes 

A(e k , uj) = — ^ 9 ; + — ; ^ ^ 9 ■ , (38) (i.e., the third-harmonic generation, e l3ut , neglected above): 



(o;-e fc ) 2 + 7 2 
2uo + ( uj 



(u; + e fc ) 2 + 7 2 
e k + ij uj 



4cj 2 + 7 2 V>-^) 2 +7 2 (^ + e fc ) 2 + 7 2 



(39) 



It is important to note the symmetries, A{c,uS) = A(— c,uj) 
and B(e,uj) = B(—e,uj), which help in the calculation of the 
total optical susceptibility. 



real 



virtual 




bleaching 



two-photon processes 



FIG. 1. Third order optical processes in neutral graphene. Bleach- 
ing: two photons are absorbed with no virtual states involved. This 
process requires emission of a photon before the second one is ab- 
sorbed. Resonant two-photon processes: a process where two pho- 
tons are simultaneously absorbed involving a virtual state. In the case 
represented here, that state is located at zero energy. These type of 
processes also occur in traditional semiconductors, when the photon 
energy is smaller than the band-gap. 27 

Analogously to P x ,i, the non-linear polarization P x ^ is ob- 
tained from 



Px,3 = iev F gsgv ^ sin [x 3 (t) - xl(t)] . 



(40) 



Replacing the expression for x$(t) in Eq. (40) we obtain (for 
neutral graphene at zero temperature) 



x,3 



(co;bl) E -icut 

Xs 2 6 



(u;;2 7 ) E 0-iut 

Xs 2 



where an d X^' 27 ^ are given by 



(oj;bl) 

Xs 



(w;2 7 ) 
X3 



34 

-aoVViTra- 



' huj 2r y 2 



-cr /3 7 



-(JoWiTra 



3vi 



2uj 



Huj 3 4uj 2 + 7 2 



h.c, (41) 
(42) 

-a Pj2, (43) 



with the respective imaginary parts being negligible in the 
regime fry <C huo. Note that both Sftx^' 6 ^ and ^Xs^ 
are negative due to saturate absorption. Also, both pro- 
cesses contribute to the imaginary part of the refraction in- 
dex of graphene. We should note that the limit 7^0 



(u;bl) 



and 



can be taken in fi u but not in /? 7 . Indeed, 5ftx 

SRxi^' 27 ^ correspond to two different physical processes typi- 
cal of semiconductors pZMl bleaching and virtual two-photon 



X 3 



(3w) , AJ Snvj? 
v -= cr Wia - 



Ahuj A 



(44) 



The transmittance of free-standing graphene for normal in- 
cidence is obtained from 



r = 



(45) 



|l + xM/(2e c)| 2 ' 
Taking into account the non-linear corrections, we have 

9?xM = <t (1 - Pi) ~ °oPu,/2 • (46) 

Since ao/(2eoc) = na, the transmittance of neutral graphene 
at zero temperature is 



T — 1 — ira + 7ra/3 7 + 7ra/3 UJ /2 , 



(47) 



where we have expanded Eq. (45 ) in the small parameter a. 
Higher-order terms are negligible except for very high field 
intensities, that is, /3 7 ,/?u; ^ 1- As expected, the non-linear 
contributions, /3 7 and (3 U , induce a higher transparency of 
graphene, which increases as Wi also increases. [The imag- 
inary part of x(cj), which we have neglected, gives a small 
correction to Eq. ( [47] ) — see the top right panel of Fig. [2] for 
the magnitude of the imaginary part of ^x(^) m m e rotating 
wave approximation.] 

It is important to realize that although the value of Wi, 7 , 
coming from /3 7 , suggests that the absorption of graphene 
would saturate for moderate intensities, the much larger value 
of Wi iU} , coming from f3 u , shows that graphene still absorbs 
light due to virtual resonant two-photon processes, even in the 
event of negligible bleaching. This is possible because, as al- 
ready noted, the two processes -bleaching and the two-photon 
process- excite different electronic states. If light of broad 
spectral range is considered, instead of monochromatic light, 
the analysis will be more complex than the one presented here. 



C. Approximate calculation of the non-linear susceptibility to 
all orders in the intensity of the field within the RWA 

In the previous section we made a perturbative calcula- 
tion of the non-linear optical susceptibility of graphene, up 
to first order in /3 7 and /3 W . The exact calculation to all orders 
is not possible. However, an approximate calculation of the 
non-linear susceptibility valid to all orders in Wi can be ob- 
tained using the rotating wave approximation (RWA). Within 
the RWA, the solution of X(t) = X is written as 



X(t) 



xe 



-iujt 



ye 



iuot 



(48) 
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Inserting X(t) in Bloch's equations, we obtain 

• i ( i • • \ -dkEo AT 
x + (72 + £e — iuj)x — 1— — N , 

V + (72 + ie + uo)y = -i^^-N , 



(49) 
(50) 



N + 7l jV « 4^9(x -y). (51) 

The explicit solution of the above set of equations is obtained 
by series resummation (assuming, for simplicity, 71 « 72 = 
7) and reads: 

r <ifc£ 7 + 2(0; — 2Virfc) 



y = iN 

N : 



2a; (u;-2v F /c) 2 +7 2 
dfc^o + 2-upfc) - 7 

2a; (u; + 2v F /c) 2 + 7 2 
no 



l + (J(fe) 
where 5(fc) is given by 



(52) 

(53) 
(54) 



<S(fc, 0) = /3 7 2 sin 2 - 



(uj - 2v F k) 2 + 7 2 



+ 



(cj + 2v F /c) 2 + 7 2 y ' 
and P is defined as fi = (8/3) /3 7 . 



(55) 



0.9 



0.6 



— co=0.3 eV, r 2 =10 meV, T=30 K 

- - co=0.3 eV, r 2 =10 meV, T=300 K 
■ — perturbative result 
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we obtain 

Pit) = ^e-^VH + ix"H] + c.c. . (57) 

Considering zero temperature, with fi > denoting the chem- 
ical potential, the susceptibility reads 



7Tp J2fi/h-f 



(58) 



■My) 



and 



7 tt(3 7 A 



y 



l + y 2 +uj 2 /j 2 



(59) 



where 



9(y) 



l 



l 



G,-^/ 7 ) 2 + i Q, + ^/ 7 ) 2 + i 



(60) 



To zero order in j3 we have the usual results (intra-band 
contributions excluded): 



X'(cj) = 1 + 



E 

s=±l 



S H(jJ 

— arctan — 
7T 7^ 



2//S 



(61) 



and 



1 (2/i + hw) 2 + ft 7 2 
" a °2^ (2/i-M 2 + ^7 2 



(62) 



To first order in f) and for neutral graphene at zero tempera- 
ture, we obtain for the real part of the susceptibility the ap- 
proximate result 



FIG. 2. Optical susceptibility of graphene as function of j3 within the 
RWA to all orders in Wi. In the left panel we plot the real part of the 
optical susceptibility of graphene as function of the dimensionless 
parameter f3 and in the right top panel we plot the imaginary part 
of the same quantity (/x = 10 -3 eV). The transmittance, T, of light 
through graphene at normal incidence is plot in the right lower panel. 
The vertical dashed line corresponds to /3 7 = 1. The results are for 
huo = 0.3 eV and fry = 10 meV and two different temperatures, 
T — 30, 300 K. The dashed-dotted line is the perturbative result 
given by Eqs. ( [46] ) and ( [47] ). 

The calculation of the polarization follows the same proce- 
dure as before, and using 



dO 



sin 



2tt vTT^ 2 - 1 
1 + a 2 sin 2 ~ a 2 VTT^ 2 ~ ' 



(56) 



(63) 



It is apparent from Eqs. ( [46] ) and ( [63] ) that the contribution 
coming from the bleaching process is exact [first term in 



Eq. (63)], whereas the contribution from virtual two-photon 
process is overestimated by a factor of two in the RWA. 

In Fig. [2] we compare the perturbative (dashed-dotted line) 
results given by Eqs. ( [46] ) and ( [47] ), left and right-bottom pan- 
els, respectively, with the RWA, valid for an arbitrary value of 
Wi (we plot results for two different temperatures: T = 30 
K, solid line, T = 300 K dashed line). Clearly, the perturba- 
tive result and the RWA calculation agree well up to f3 1 = 1, 
which sets perturbation theory validity limit. Beyond that 
value a non-perturbative approach is necessary and one has 
to rely on the RWA for drawing quantitative conclusions. 
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IV. AN EFFICIENT GRAPHENE-BASED 
PHOTODETECTOR 

In this section we describe the absorption of light by a de- 
vice composed of optical cavities and a single graphene sheet. 
Only the linear optical susceptibility will be considered, ex- 
cept when the light intensity inside the cavity is of the order 
of Wi, 7 (we remark that for telecommunication devices and 
photodetectors this will hardly be the case). 



A. Properties of a mirror and of an empty optical cavity 

We start with the well-known case of an empty optical cav- 
ity. This allow us to introduce important concepts and fix the 
notation. The scattering matrix of a partially-silvered mirror 
is characterized by two pairs of reflectance and transmittance 
coefficients; such a mirror is shown in Fig. [3] 



E, 



Ei 



I 



J?ij{ _|_ graphene 



(t a ,r a ) 



I 



E o,+ 
>- 



Eo r 



FIG. 3. Graphene inside an optical cavity. The cavity is defined by 
two equal mirrors with reflectance amplitudes r a and r&, and trans- 
mittance amplitudes t a and tb. 

The S— matrix relates the amplitudes of the incoming 
waves, Ei i+ and Ea,-, to the amplitudes of the outgoing 
waves, Ei - and Ea,+, according to 



(64) 



On the other hand, the transfer matrix M m relates the fields 
on the two sides of the mirror according to 







t a n 




E i,+ 






_r a t b 







E A,+ 


1 







t a h ~ r a r b r b 
-r a 1 



Ei,+ 



(65) 



Let us now discuss few properties obeyed by the reflectance 
and transmittance amplitudes. 30 From the conservation of the 
energy flux 



\Ei,+ \ 



\Ea,- 



\Ea,- 



we find 



From t a r* b 



\ta\ 2 + \Va\ 

\h\ 2 + \n\ 

t a r* b = 
-t b r a we find 



= 1. 



(66) 

(67) 
(68) 
(69) 

(70) 



which can be used to show that \t a \ = \tb\ and \r a \ = \r b \. 
Furthermore, the determinant of the transfer matrix reads 
detM m = t a /t b , implying that |detM m | = \t a /t b \ = 1. For 
systems with inversion symmetry, we can write t a = t b = t 
and r a — r b = r. Thus, the S— matrix reads 



S 



t r 
r t 



(71) 



and the transfer matrix can written as 



1 


t 2 — r 2 r 




"t + |r| 2 /r 


t 


-r 1 




-r/t 



1/t* r/t 
r*/t* 1/t 



(72) 



Writing r = \r\e tar and t = \t\e iat , the relation tr* = —t*r 
implies that e^^^-at) _ is, a r = a t ± tt/2. Using 

these last relations, M m can be written as 



-1 



1 



(73) 



If two of these mirrors are separated by a distance L we 
have to define the transfer matrix associated with the free 
propagation from the first to the second mirror. Since 



E+(x,t) =E + e i ^ x ~ UJt \ 
E_(x,t) = E_e~^ kxJrUJt) , 



(74) 
(75) 



then at a distance L to the right the E+ (x + L,i) has an extra 
phase of e tkL whereas the E_ (x + L, t) has an extra phase of 
e~ lkL . Thus we have 



or 



(76) 







- e ikL Q 






Ea-,l 




e~ ikL 







^4,+ 




" e -ikL Q 




Ea,+,l 


_ E *,- _ 




e ikL 




Ea-,l _ 



(77) 



where Ea,+/-,l represents the amplitude of the for- 
ward/backward propagating field at the right end of the cavity 
and the matrix 



M f (L) 



-ikL 









ikL 



(78) 



defines the free propagation to the right. Then, the transmit- 
ted field through two mirrors at a distance L from each other 
follows from 



Ei. 
E L 



Explicitly, we have 



Ei,+ = 



M m ■ M f (L) ■ M„ 



Eo,+ 
Eo,- 



-ikL 



[1 - (r*Ye MkL ]E ,+ , 



(79) 



(80) 
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writing r* = \r\ 



we obtain 



T - 



\Eo, 



l^, + l 2 



|£| 4 + 4|r| 2 sin 2 (/cL) 



(81) 



since a r = it, as implied by Fresnel equations. Thus we 
have perfect transmission for kL = nir or A = 2L/n, with 

n = 1,2, 3, The longest wavelength for which perfect 

transmission is possible is A = 2L. It is straightforward to 
show that the transmission is strongly suppressed for other 
choices of kL. For instance, when kL is a multiple of tt/2, 
we have 



T 



\t\ 4 



t 



•4|rl : 



4|rl : 



< 1. 



(82) 



In the above, we have admitted high-quality mirrors, \r\ ^> \t\, 
to simplify the denominator. 

The introduction of a graphene sheet inside the cavity leads 



to light absorption and the relation (81 ) is modified. In the 
following section, we demonstrate how to explore the physics 
of an optical cavity to devise an efficient graphene-based pho- 
todetector. 



B. Graphene in an optical cavity 

We describe the transmission of light through a graphene 
sheet inside an optical cavity taking into account the linear 
optical- susceptibility of graphene. 



transmittance 



50 
40 

CD 

I 30 

CD 
O 

5.20 
10 



rr 



T3 



i .I i I i j u I i l. i 





0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 
kl_/7t kl_/7i kU% 




0.9 0.95 1 1.05 1.1 0.9 0.95 1 1.05 1.1 0.9 0.95 1 1.05 1.1 
kl_/7t kl_/7l MJ% 



where 2rj = Z xM, with rj = rj' + if]", and Z ~ 376.7 Vt 
is the vacuum impedance. For neutral graphene at zero tem- 
perature r] is essentially a real number for frequencies below 
the visible spectral range. The transmission through the cavity 
with graphene at position x g and the second mirror at position 
x = L follows from 



Ei 



,+ 



= M m -M f (x g )-M g -M f (L-x g )-M n 



The matrix 



M = M m • Mf(x g ) • M g • Mf(L 



• Mrr 



(84) 



(85) 



is the full transfer matrix of the device. The transmittance and 
the reflectance are defined as 



T 



K 



|Mn| 2 



M 21 



M 



li 



(86) 
(87) 



respectively, and M\\ and Mi 2 denote the matrix elements of 
M. We note that 1Z ^ 1 — T due to absorption by graphene 
(we are assuming lossless mirrors). The absorbance is defined 
'&$A=l-K-T. 
For x g = L/2 we have 

T = 

K = 



t 



■ 77 — 2r]\r\e lkL 

j_ ^)| r | _ /y(l -f 



(1 



^^| r |2g2i/eL|2 



(88) 



\r\ 2 )e ikL 



(1 — 77) |r*|e 



2ikL\2 



\1 + T] - 2rj\r\e ikL - (1 



, jj\ | r |2g2i/cL|2 



489) 



In the limit r — > we recover the well known result, T 



1/| 1 + ??| 2 [see also Eq. d45b] and K 



\v\ 2 /\i + v 



2[24] 



The effect of graphene in the cavity is to reduce the inten- 
sity of the odd-orders (n = 1, 3, 5, . . .) of perfect transmis- 
sion in the otherwise perfect cavity. From Fig. [4] it is clear 
that the reduction of transmission of the odd orders is divided 
between reflection and absorption, the latter taking the major- 
ity of the incoming power. The transmission is still unity for 
even-orders (n = 2, 4, 6, . . .). In Fig. [5] we show the depen- 
dence of the absorbed power as function of the transmittance 
t 2 of a mirror. Clearly, this dependence is not monotonous, 
displaying a maximum around t 2 ~ 0.045. An analytical ex- 
pression for the value of t 2 for which the absorption is maxi- 
mum can be readily obtained from Eq. i 



FIG. 4. Transmission spectrum for a Fabry-Perot cavity with 
graphene L/2. The mirror transmittance is t 2 = 0.045. If we tale 
A = 1000 nm, then the range kL/n in the lower panels spans the 
wavelengths from A = 900 nm to A = 1100 nm. We have taken 
fir/ = 7 meV and T = 300 K. Regarding the choice of A, we note 
that the first HeNe laser was working at the spectral wavelength of 
1 150 nm, that is in the infrared. 

We write the transfer matrix of graphene as 22 



C. Double optical cavity 

The goal of the present section is to discuss a device able to 
enhance light absorption relatively to the single cavity system 
discussed above. To that end we consider a graphene sheet 
inside an optical cavity of length L (as in the previous sec- 
tion) followed by an empty quarter- wavelength cavity, as rep- 
resented in Fig. [6] 

Placing graphene at the middle of the first cavity and choos- 
ing its length L such that L = A/2, where A is the wavelength 
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100 
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transmittance 
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FIG. 5. Absorbance of graphene for the single cavity system (see 
Fig. [3}. Left: Absorbance as function of mirror transmittance 
t 2 G [0,0.5]. Top right: Absorbance of graphene as function of 
reflectance. Bottom right: reflectance of graphene as function of 
transmittance. The calculations assumed T = 300 K. 



FIG. 7. Absorbance of graphene for the double cavity setup (see 
Fig.(6jl. Left: Absorbance as function of mirror transmittance t 2 £ 
[0,1]. Top right: Absorbance as function of reflectance. Bottom 
right: Reflectance as function of transmittance. The calculations as- 
sumed T 300 K. 



kL = 7T k L 2 = rep 



E, 



(t a ,r a ) 



Ej± + graphene 



Ea- 

(t b ,r b ) 



Eq. 



FIG. 6. The double cavity system: Light impinges from the left- 
hand side. The second cavity has half the length of the first cavity, 
L 2 L/2. 



of the light, we expect that graphene will present an enhanced 
absorption at this wavelength, at least for a cavity with a high 
finesse. This intuitive picture is developed from considering, 
as a rough approximation, the formation of standing waves 
within the cavity having their maximum amplitudes at the cen- 
ter of the cavity. This picture in confirmed by simulations, as 
shown in Fig. [4] Also from Fig. [5] it is clear that there is an 
optimal value of t 2 for which the absorption can be as high 
as 50% (the quantitative results are robust for small changes 
graphene 's position relatively to the center of the cavity). 

The absorbed intensity can be pushed up to ~ 100% by 
building an optical cavity containing graphene, followed by an 
empty quarter- wavelength cavity, that is, a cavity with length 
L 2 such that L 2 = A/4. This setup leads to an enhancement 
of the absorption, which is about twice as large as that of the 
single cavity setup, as can be seen in the left panel of Fig. [7] 
The dependence of the absorption on the wavelength is shown 
in Fig. [8] 

A physical qualitative argument for the absorption enhance- 
ment effect in a double cavity is reminiscent of a quantum 
particle in a box with a permeable wall. Let us consider first 



a box with origin at x = and length L + L/2. A wall is 
located at x — L. If the wall is impermeable, the fundamental 
mode of the the first box is A = 2L and that of the second box 
is A = L. If the wall becomes permeable the two modes hy- 
bridize and in the ground state the probability density grows 
in the first box at the expenses of the probability density in 
the second box. Translating this into our problem, the quarter 
wavelength cavity interference between the wave reflected by 
the third mirror and the forward propagating wave effectively 
suppresses the transmission at wavelengths A = 2L, forcing 
the photon to spend more time in the first cavity, and thus in- 
creasing the absorption by the graphene sheet. 

It is worth stressing that the maximum of absorption takes 
place for cavities with \t\ 2 ~ 0.1, a convenient figure from the 
point of view of micro-fabrication since not much effort has 
to be put on building highly reflective mirrors. 

From a theoretical point of view, the calculation of the prop- 
erties of the two coupled cavities follows from the transfer ma- 
trix method, reviewed in |IVB| As in Eq. (84), the incoming 
and outgoing field amplitudes are related as 



Ei,+ 
E, 



Eo,+ 
E o,- 



(90) 



where the transfer matrix of the two cavities is given by 

M cav = M m -Mf{x g )-M g -M f {L-x g )-M m -Mf{L 2 )-M m . 

(91) 

In Figs.|7]and[8] we have considered x g = L/2. In this case the 
field amplitudes for A = (2n + l)L/2 (with n = 0, 1, 2, . . .) 
have a maximum at the center of the cavity. As for the case of 
the single cavity, it is possible to derive analytical expressions 
for both T and 1Z; we obtain 



T 



w 



|1 + 7? + Aie 



ikL 



A e 3ifcL|2 



(92) 
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FIG 8. Absorbance as function of kL/n = 2L/X for the double 
cavity system (t 2 = 0.09). If we consider A = 1000 nm, then 
the range of kL/ir in the lower panels spans the wavelengths from 
A 900 nm to A = 1100 nm. 

In the above, A 1 = -[77(2 + |r|) + \r\]\r\, A 2 = [77(1 + \r\ + 
|r| 2 ) - \r\]\r\ and A 3 = (1 - rj)\r\ 2 , and 

n = ^ |(1 + 77)|r| + Aie ikL + A 2 e 2 ^ L + A 3 e 3 ^ L | 2 , 

(93) 

withAi = -[77(l + |r| 2 + |r| 3 ) + |r| 3 ], A 2 = (^-1+2^)^1 
and A3 = (1 — rj)\r\. As before, taking the limit r Owe 
recover the well known values for T and 1Z in the absence of 
mirrors. 

It is interesting to note the hysteresis in the absorbance 
as function of the reflected power — some values of the re- 
flectance admit three possible absorbances; see top right panel 
of Fig. [7] Each point on the A versus 1Z curve corresponds to 
a given value of the mirror transmittance, |t| 2 , which therefore 
can be viewed as an external driving field. 

Finally, we note that about 100% absorption can also be ob- 
tained for a single cavity if the second mirror reflects 100% of 
the impinging light. The curves A, TZ, and T as function of t 2 
are different in this case, however, from those given above; in 
particular, A shows no hysteresis (see note after Conclusions 
section) 



V. CONCLUSIONS 

In the present article, we have developed the theory a 
graphene-based photodetector with nearly 100% efficiency for 
photon frequencies around a predefined value. The proposed 
setup is general and should work in a vast spectral range. 

We have also clarified the role of the non-linear optical 
susceptibility in determining the properties of the cavity- 
graphene system. We have shown that non-linear terms are 
irrelevant for moderate light intensities. 

The most efficient photodetector is built from combining 
a half- wavelength cavity (size L = A/2) followed by a sec- 
ond quarter- wave length cavity (size L 2 = A/4). This system 
improves the absorption by a factor of two relatively to the 
single cavity. As noted above, the two-cavities photodetector 
has about the same absorbance as a single cavity with the sec- 
ond mirror having zero transmittance and the first one having 
t 2 ~ 0.9. 

If real-time control of the mirrors transmittance is feasible, 
then it will be possible to obtain hysteresis in the absorbance 
for the same reflectance value. Whether this can be used as an 
optical logical gate is so far unclear and will be left for future 
research. 

Although in the schematic figures for the cavities graphene 
appears floating in the air, in practical terms it will be de- 
posited on a dielectric. The mirrors can also be made of di- 
electric materials, the so called Bragg mirrors. There are com- 
puter codes for simulating mirrors with the prescribed optimal 
value of t 2 given in the text. The setup itself can be built by 
micro-fabrication using standard techniques. 

Note: during the final state of writing, we became aware 
of an experimental paper, entitled ''Microcavity -integrated 
graphene photodetector". In that work a single cavity de- 
tector has been built. 31 Our theory is a full analytical account 
of the physics of two similar devices, one of them having the 
same theoretical efficiency as the device of that paper but a 
different qualitative response as function of the amplitude t 2 . 
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